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Investigation of transfer processes in dense weakly or fully ionized many-component 
gases is topical in view of construction and improvement of new gaseous lasers, plasmo- 
chemical reactors for ozone synthesis, and air cleaning of both nitrogen and carbon oxides 
and chlorine compounds. 

A study and calculation of transfer coefficients taking into account the interparticle 
interaction nature is one of the major problems in this investigation. Recent papers by 
Murphy [1-4] are devoted to these calculations for some special low-density mixtures bas- 
ing on the Boltzmann kinetic equation for point particles and using accurate interparticle 
potentials of interaction. Nevertheless, there are some problems in calculation of trans- 
fer coefficients for high-density mixtures. They are caused mainly by the fact that the 
Boltzmann kinetic equation is suitable for rarefied gases and plasmas only. 

In present work a calculation of transfer coefficients for many-component dense gases 
for charged and non-charged hard spheres is carried out using the Enskog-Landau kinetic 
equation which takes into account realistic particle sizes. 

The Enskog-Landau kinetic equation was obtained in [5,6] for a one-component system 
of charged hard spheres in the electron compensation field starting from the Bogolubov 
hierarchy with modified boundary conditions in the pair collision approximation. In paper 
[7] this equation was generalized for the case of many-component system and usually reads: 
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where f a {x a \t) is the one-particle distribution function of a-component (x = {r,p}), M 
is a number of components. 

The collision integral in the right-hand part of (1) consists of several terms. Such a 
structure is produced due to the additivity of the interaction potential which is modelled 
for the mixture of charged hard spheres as a sum of the short-range part (hard spheres) and 
the long-range one (Coulomb particles). 1^ and 1^ are the zeroth and first expansion 
terms of the Enskog collision integral with respect to the spatial inhomogeneity [5]. Imp is 
the collision integral of the kinetic mean field theory KMFT [8,9]. It is the first order on 
the long-range interaction. The last term 1^ is the generalized Landau collision integral 
and it is the second order on the long-range interaction. 

The kinetic equation (1) can be solved with operating the standard Chapman- Enskog 



method [10]. The unknown distribution function in the first approximation can be taken 

as 
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where is the zeroth approximation. Usually, it is chosen as the local-equilibrium 
Maxwell distribution function. <p a is a correction which is expressed through the Sonine- 
Laguerre polynomials [10]. 

Having the solution to the kinetic equation (1), one can calculate the stress tensor 
and the heat flow vector in a system. These quantities are expressed via such transport 
coefficients as bulk k and shear rj viscosities and thermal conductivity A, respectively: 
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where a ab = (a a + u b )/2, a a , a b - hard sphere diameters, glf - quasiequilibrium binary 
correlation function, n a - particle number density, ji ab = m a m b / (m a + m b ) - reduced mass, 
A a (0), A a (l) and B a (0) are zeroth and first coefficients of expansion of the correction a 
on the Sonine-Laguerre polynomials. 

In numerical calculations of transfer coefficients (3)-(5) for some binary and ternary 
mixtures of neutral and charged hard spheres we considered their temperature and concen- 
tration ratio dependences. The values of hard sphere diameters were fixed: <7h c = 2.15A, 
o"Ar = 3.405A, o"Kr = 3.67A, axe = 3.924A. These values were borrowed from [11]. 

Calculating A a (0), A a (l) and B a (0) we faced with the problem of divergency of the 
so-called f2-integrals at large distances. To avoid this circumstance we should change 
upper limit of integration and use finite screening radius D instead of infinity. This 
screening radius, in contrast to the Debye formula for point-like particles, explicitly takes 
into account realistic sizes of charged particles in accordance with [12]. Contact values of 
the binary correlation function in expressions (3)-(5) were borrowed from [13]. 

Not only the transport coefficients k, rj, and A were studied. Numerical calculation for 
thermal diffusion D% and mutual diffusion coefficients D a/3 has been performed as well. 
For two- and three-component mixtures of neutral and charged hard spheres we stud- 
ied their dependences on density, temperature, and concentration ratio of some mixture 
components [14]. Figs. 1 and 2 show our calculations in brief. 

The obtained Enskog-Landau kinetic equation for charged hard spheres turned out 
to be very useful for several purposes. First of all, the collision integral of this equation 
does not contain a divergency at small distances. Secondly, the normal solution and all 
transport coefficients have analytical structure. They can be easily used to study some 
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Fig. 1: a) - temperature dependence of shear viscosity 77 for mixture 40%H2, 60% Ar at pres- 
sures P = 300, 400, 500 atm., respectively. Curves are plotted in order of increasing pressure. 
Though this model is not entirely suitable as a hard spheres model, because the H2-molecule 
is not spherical, nonetheless obtained results give an information about the successfulness of 
the application of the model for nonspherical molecules. Somewhat worse agreement is between 
the theory and experimental data [15] when temperature increases. At high temperatures the 
role of the attractive part of interparticle interaction potential becomes weaken in comparison 
with the repulsive one. However, this behaviour can be explained by ignoring in the theory the 
temperature and concentration dependence of hard sphere diameters for molecules of Ar and 
H2. At the same time it is known that a depends both on T and n. b), c) and d) - temperature 
and concentration ratio dependences of shear viscosity r/ and thermal conductivity A coefficients 
for the mixture Ar + -Kr + at the fixed total density n = 2 ■ 10 26 m~ 3 . These calculations were 
performed for three fixed concentration ratios of the heavy component x Kr + = 0.25; 0.5; 0.75, 
iKr+ = n Kr+ 1 ( n Ar+ + n Kr+ ) • The coefficients r/ and A have very similar temperature depen- 
dence. But concentration ratio dependence for A has a different curve fit. More detailed analysis 
is given in [14]. 
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Fig. 2: Temperature dependences of transport coefficients: shear viscosity r] and thermal 
conductivity A a), and mutual diffusion coefficient D^b and thermal diffusion coefficient 
D T b) for the system 26.7%Ar+, 73.3%He at n = 1 .0487 • fO 27 m" 3 , or A = 0.00976. 

specific systems. Finally, the analytical structure of transport coefficients allows us to find 
fast and easily systems, which can be best described by the obtained kinetic equation, 
as well as density and temperature ranges, where the agreement between the theory and 
experimental data is the closest. 

A very important step in this theory is to calculate a dynamical screening radius in a 
system. Partially this problem has been already solved in our recent paper [16]. 
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